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Abstract. We present the results of our study of the variabihty time scales of a sample of 27 long Swift Gamma-Ray Bursts 
(GRBs) with known redshifts. The variability time scale can help our understanding of fundamental ORB parameters such as 
the initial bulk Lorentz factor and the characteristic size associated with the emission region. Fast Fourier Transform (FFT) 
techniques were used to extract a noise threshold crossing frequency, which we associate with a variability time scale. The 
threshold frequency appears to show a correlation with the peak isotropic luminosity of GRBs. 
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INTRODUCTION 

The time variability in Gamma-ray Bursts (GRBs) is crucial to our understanding of the characteristic size associated 
with GRBs. 

In 2000, Fenimore and Ramirez-Ruiz first proposed a correlation between variability of GRBs and peak isotropic 
luminosity [1]. Since then a number of authors have provided further support for this correlation 13, HI S, S S IS]- 
However, these authors used a variety of definitions for variability with various parameters and smoothing methods. 
The lack of a universally accepted definition for variability is a major short coming and poses problems in comparing 
and evaluating the results of previous studies. 

In our analysis we have used Fourier analysis techniques to probe various frequencies and their strength in 
GRB light curves. In this paper we associate a threshold frequency, which is the frequency where the signal crosses 
the noise level, as a potential variability indicator. 

Due its fast slewing capability, the Swift Gamma-Ray Burst mission [8] has enabled more redshift measure- 
ments of GRBs than before. With its highly sensitive primary instrument, the Burst Alert Telescope (BAT) [9], Swift 
provides high quality GRB data on which a Fourier analysis can be performed and a variability time scale extracted. 
Our sample contains 27 Swift BAT GRBs with good spectral measurements. We extracted the threshold frequency 
for all the GRBs and show this parameter is correlated with the isotropic peak luminosity. 

METHODOLOGY 

The discrete Fourier transform of a series of N counts or numbers is given by 
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To calculate the power spectrum, the Fast Fourier Transform (FFT) algorithm was used OlOll . The power spectrum is 
defined as [11] 
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We have used BAT event-by-event data and utilized the "IDL Extract" softwar^H package to generate power spectra 
for our analysis. We divided the burst into 2 segments and averaged the corresponding power spectra to get the 
final spectra. A typical power spectrum from a GRB is shown in Fig. [T] The low frequency power-law component, 
sometimes referred to as "red noise", represents the signal from the source. The high frequency region, called "white 
noise", is background. The threshold frequency, /th, is the frequency at which the red noise crosses the white noise 
level. In order to extract the threshold frequency, the power spectra were fitted by a broken power law 
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FIGURE 1. The power as a function of frequency. A broken power-law is fitted to the power spectrum to obtain the threshold 
frequency where the red noise (source signal) intersects the white noise (background). 

When calculating the isotropic peak luminosity (Liso), we need to take into account the fact that in the rest frame of 
the GRB, photon energies are higher than those in the observer frame. The observed peak flux for the source-frame 
energy range l.OkeV to 10,000keV was calculated using observed spectral-fit parameters for the GRB sample and 
appropriately z-corrected energy limits in the integral. The isotropic peak luminosity is related to the calculated flux 
as follows: 



Liso =47r^/^Fobs• 
Here Fobs is the observed peak flux and the luminosity distance, d^, is given by. 
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For the cuiTent universe we have assumed, = 0.27, £2^ = 0.73, Q.^ = 0.0 and a Hubble constant (Hq) of 
70kms^'Mpc^' = 2.268 x 10^'**s^'. The redshift measurements were taken from online archives of the Gamma- 
Ray Burst Online Index (GRBO?© and verified by using the GCN circulars^ 



' http://idlastro.gsfc.nasa.gov/ftp/contrib/rxte/ 
" http://lyra.berkeIey.edu/grbox/grbox.php 

^ Gamma-ray burst Coordination Network jhttp://gcn.gsfc.nasa.gov) 



RESULTS 



From a sample of 100 Swift BAT GRBs with spectroscopically confirmed redshifts, we selected 27 GRBs with good 
spectral measurements. We analyzed event-by-event data of this sample and obtained power spectra and fitted these 
with a broken power law and extracted the threshold frequencies. We find the average a for the sample is 1 . 12 ± 0.05 
and j3 is consistent with zero for all fits. In Fig. |2] panel (a), we show the calculated isotropic peak luminosity 
as a function of the extracted threshold frequency (with the appropriate z correction). As seen in the figure, the 
threshold frequency appears to be correlated with the isotropic peak luminosity. The Pearson's correlation coefficient 
is 0.69 ± 0.03, where the uncertainty was obtained through a Monte Carlo simulation. The probability that the above 
correlation occurs due to random chance is 3.4 x 10^^. Our best-fit yields the following relation between Liso and /th: 

logLiso = (52.0 ± 0.2) + (1 .4 ± 0.2) log(/th (z + 1 )). (5) 

Panel (b) in Fig.|2]shows a histogram of the threshold frequencies. The lowest redshift-corrected threshold frequency 
in the sample is approximately 0.2 Hz and the largest is about 20 Hz. Our results imply the smallest variability time 
scale to be approximately 50 msec. 




DISCUSSION 

It is important to investigate how the observed brightness of GRBs affects the extracted threshold frequencies because 
the threshold frequency also appears to be correlated with the observed brightness (see Fig. |3] panel (a)) and this 
potential observational bias needs further investigation. 

FFT power is roughly proportional to the brightness of the burst. Hence when a burst of a given luminosity is 
closer and therefore brighter, the red noise is enhanced above the constant Poisson noise level (white noise) and results 
in a larger threshold frequency. Given a burst of luminosity Liso and frequency /th, the same burst with a different 
luminosity LJ^q would fall on the following L!^^^ — line: 

logLLo = log(Liso//£) + alog.4. (6) 

The mean value of a for our sample is 1.12 ±0.05, which does not fully explain the slope of 1.4 ±0.2 (see Fig.|2|i 
obtained from the Liso-/th correlation. 



We note also that the uneven distribution of redshifts (see Fig |3] panel (b)) might play a role in the correlation 
with flux. Roughly half of the sample lies in the redshift range 1.5 to 3.5 where the luminosity distance changes only 
by factor of 2.8, hence, the correlation between Li^o and /th partially trickles down to the flux as well. 
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FIGURE 3. Panel (a): Isotropic peak luminosity as a function of observed flux. Panel (b): The histogram of redshifts in the 
sample. 

In conclusion, we argue that the FFT threshold frequency is a parameter that provides a suggestive measure of 
gamma-ray bursts variability. For the sample of 27 GRBs analyzed, our results imply the smallest variability time 
scale is approximately 50 milliseconds. The apparent correlation between the isotropic peak luminosity and the 
threshold frequency needs further investigation for observational bias. If confirmed, this correlation is potentially 
useful as a probe of GRB microphysics. It may also prove useful as a redshift estimator. 
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